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Abstract 

We study a simple model for a neuron function in a collective brain system. The 
neural network is composed of uncorrelated random scale-free network for elimi- 
nating the degree correlation of dynamical processes. The interaction of neurons is 
supposed to be isotropic and idealized. This neuron dynamics is similar to biological 
evolution in extremal dynamics with isotropic locally interaction but has different 
time scale. The evolution of neuron spike takes place according to punctuated pat- 
terns similar to the avalanche dynamics. We find that the evolutionary dynamics 
of this neuron function exhibit self-organized criticality which shows power-law be- 
havior of the avalanche sizes. For a given network, the avalanche dynamic behavior 
is not changed with different degree exponents of networks, 7 > 2.4 and refractory 
periods correspondent to the memory effect, T r . In addition, the avalanche size dis- 
tributions exhibit the power-law behavior in a single scaling region in contrast to 
other networks. However, the return time distributions displaying spatiotemporal 
complexity have three characteristic time scaling regimes. 
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1 Introduction 



A human being is most complex organism in the nature. A human brain 
is extremely complicate in the organs of human. That is the reason that it 
is exceedingly difficult to understand brain through individual and classical 
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mechanism. The nervous tissue of the human brain contains many billions of 
neurons, large complex cells that conduct nerve impulses from one part of the 
body to another part. Scientists have investigated neural network in the brain 
through various methods [1,2,3,4,5]. Recently, neural network was reported to 
be like the scale-free network, which is characterized by displaying power law 
distribution of the degree [6]. An well known scale- free network is Barabai- 
Albert(BA) network characterized by evolving and preferential attachment 
[7,8]. Goh et al. introduce a static scale-free network with no degree-degree 
correlation [9] . The firing pattern of the neuron is similar to the avalanche pat- 
tern of the self-organized criticality(SOC) [10,11,12,13]. The propagative size 
of neuron functioning exhibits the power-law distribution. Thus, the punctu- 
ated pattern and power-law behavior occur without fine-tuning parameter. L. 
da Silva et al. offered a simple model for brain functioning similar to the Bak- 
Sneppen(BS) model [11] with the memory effect on the lattice [14]. The same 
model is investigated on small- world network [15] by Lin and Chen [16]. We 
consider an avalanche dynamics of an idealized neuron function on a uncorre- 
lated random scale-free network. The structure and the dynamics of neurons 
describe as follows. The anatomical unit of the nervous system is the neuron. 
The brain possesses about 10 10 — 10 12 neurons. Neurons are quite complex, but 
each of these is made up of dendrite, a cell body, and an axon. A dendrite con- 
ducts signals toward the cell body. The cell body is the part of a neuron that 
contains the nucleus and other organelles. An axon conducts nerve impulses 
away from the cell body. Dendrites and axons collectively are called neuron 
fibers. A nerve impulse is the way a neuron transmits information. When an 
axon is not conducting a nerve impulse, the resting potential indicates that 
the inside of an axon is negative compared to the outside. In contrast, if an 
axon is conducting a nerve impulse, an action potential (i.e., electrochemical 
change) travels along a neuron. As an axon is stimulated by an electric shock, 
threshold may be reached for an action potential. A fiber can conduct a vol- 
ley of nerve impulses because only a small number of ions are exchanged with 
each impulse. As soon as an impulse has passed by each successive portion of a 
fiber, it undergoes a refractory period during which it is unable to conduct an 
impulse. This ensures a one-way direction of the impulse. During a refractory 
period, the sodium gate cannot yet open. Every axon branches into 10 3 — 10 4 
fine terminal branches called a synaptic bulb . Each bulb lies very close to the 
dendrite (or the cell body) of another neuron. This region of close proximity 
is called a synapse. At a synapse, the membrane of the first neuron is called 
the presynaptic membrane, and the membrane of the next neuron is called the 
postsynaptic membrane. Transmission of the nerve impulse from one neuron 
to another takes place across a synapse. In humans, synaptic vesicles release a 
chemical, known as a neurotransmitter to receptors in the postsynaptic mem- 
brane increases the chance of a nerve impulse (stimulation) or decreases the 
chance of a nerve impulse (inhibition) in the next neuron [17]. The brain as a 
whole is a system capable of auto- regulations [14]. 
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Fig. 1. The log-log plot of the probability distribution function P(s) of the i?o(7) 
avalanche size as a function of the avalanche size s at the critical fitness on UCM 
for 7 = 2.1,2.5,2.9,3.0,4.0, and 5.0 with N = 10000, m = 3, and T r = 1. The 
critical threshold barriers where B = 0.039,0.080,0.121,0.130,0.212, and 0.260 
respectively. 



Model and Simulation Method 



To make a modeling of this neuron dynamics, we make networks. The networks 
is composed of uncorrelated random scale-free network. Here, the uncorrelated 
random scale-free network is called by the uncorrelated configuration model 
(UCM) [18]. We generate iV-cells by the static method. Let's make iV^-cells of 
the degree k satisfying the degree distribution P(k) ~ fc -7 . Select two nodes 
randomly and connect them if they are not connected before. We exclude du- 
plicated connections and self-connection. The network generated by the UCM 
is not only fully connected and have but also no degree-degree correlation 
Dnn{k), defined as the average degree of the nearest neighbors(NN) of the 
nodes with degree k and clustering correlation C(k), mean by the probability 
that a node of degree k form loop with two NN [18]. The number of minimal 
degree is fixed at m=3 to prevent alteration of dynamics from dangling node. 
Now, we introduce the evolution rule of neuron dynamics. Each node of the 
UCM represents a neuron and a link between two nodes represents a synapse. 
The uniform random numbers are distributed over each node. This random 
number is called by a barrier, which is the possibility of the firing of a neuron's 
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Fig. 2. Dependence of the avalanche size exponents as a function of network degree 
exponents 7 with N = 10000, m = 3, and T r = 1. 

spike. The lower barrier has a higher potentialities to fire [14] . To simplify the 
firing is occurred at an neuron with lowest barrier sequentially. And directly 
connected neurons also is fired because they are enough stimulated to fire. 
The selected neuron with the lowest barrier can not fire if the elapsed time 
after the firing is less than the refractory period T r . In fact, neuron's spike 
transfer directly through their synapse. However, the evolution of the signal 
is just considered to spread out to all of nearest neighbors. If this process is 
iterated the system reaches to critical stationary state, which all the barriers 
are above the B c barrier so-called a critical threshold. That is, neural sys- 
tem is self organized without well tuning parameter to stationary state. And 
the brain functioning is occurred abruptly underlying avalanche dynamics. 
The purpose of this research is to investigate how the avalanche dynamics is 
changed according to varying degree exponent 7 of networks with different re- 
fractory periods T r . For investigating spatiotemporal correlation of avalanche 
dynamics, we examined first and all return time distribution(RTDs). 



3 Results 



In stationary state, we consider the branching process of a avalanche dynamics 
to be unrestricted from network size corrections [19]. The avalanche is always 
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Fig. 3. The log-log plot of the probability distribution function P(s) of the Bo(T r ) 
avalanche size as a function of the avalanche size s at the critical fitness on UCM 
from T r = 1 to T r = 10 with N = 10000, m = 3, and 7 = 3. The critical threshold 
barrier is Bq = 0.13 for all T r . 

started from the hub neuron and all of the neurons with Bi > B , where B is 
an auxiliary parameter, are treated as inactive neurons. A B avalanche size S 
is defined as the number of the firing less than B . As B —>■ B c , the avalanche 
size distribution follow a power-law behavior P(S) ~ S~ T with an exponential 
cutoff. Figure 1 shows avalanche size distribution for the -80(7) avalanche 
with the degree exponents of the networks, 7 = 2.1,2.5,2.9,3.0,4.0, and 5.0 
The avalanche size distribution follows power-law behavior, P(S) ~ S~ T ^ 
extending over more larger regime than on other scale-free networks such 
the BA network [7] and the static model introduced by Goh et al. [9,20]. 
The more interesting thing is that the power law behaviors of the avalanche 
size distribution do not exhibit the crossover between two different scaling 
regimes. The avalanche size distribution shows a short intermediate regime 
and follows exponential decay at the cut-off regime. The absence of the two 
regimes in the avalanche size distribution may be able to explain by following 
two reasons. One of the reasons is that an average degree (k) is not fixed with 
different degree exponent 7 as compared with another static models by the 
Goh's algorithm [9] and in addition is increased as 7 gets smaller. Another 
reason is the absence of clustering correlations C(k) as well as degree-degree 
correlations D nn {k) [18]. Figure 2 presents the basic critical exponents r, the 
so called avalanche size exponent from different 7 , where the increment is 0.1 
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Fig. 4. Dependence of the avalanche size exponents as a function of refractory 
periods T r with N = 10000, m = 3, and 7 = 3. 

in 2 < 7 < 3 and 0.5 in 7 > 3, with T r = 1. As we can observe, the critical 
avalanche size exponent is the same as mean field result i.e, r ~ 1.5 for 7 > 2.4. 
It is difficult to compare the avalanche size critical exponent r in the UCM with 
the BA network because the avalanche size distribution on the BA network 
shows different power-law behavior with two regimes [21,22]. Occasionally, the 
critical thresholds are very close value each other, f c = 0.086±l(on UCM) and 
f c = 0.089 ±2(on BA) within error bar. Even though it is not universal value, 
f c is seem to a criterion distinguishing different avalanche dynamics of 7 < 2.4 
from mean field result for 7 > 2.4. In the lattice [14] and small-world network 
[16] with low rewiring probability <fi = 0.01, the levy-flight exponents and the 
avalanche size exponents increase according to the increment of the refractory 
period T r , respectively. Figure 3 illustrates the probability distribution for 
the avalanche size on UCM with different T r . At the same time, in Fig. 4, we 
demonstrate the dependence of the exponent r from 1 to 10 with Bq = 0.13. 
All r are not changed according to varying T r unlike the results on lattice 
[14] or small- world network [16]. In the lattice, the avalanche dynamics of the 
firing is propagated further far away from first update neuron because the 
firing is rejected during S < T r as refractory period T r increase. But in scale- 
free network, the active neuron returns soon to the hub and the direct linked 
node in stationary state because a hub neuron has many nearest neighbor's 
neurons. Accordingly, the firing neuron is not evolved far from hub quickly 
with refractory period. For that reason, the memory effects of refractory time 
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Fig. 5. The log-log plot of the probability distribution function P(s) of the first 
return time as a function of the time t at the critical fitness on UCM with N = 10000, 
m = 3 at T r = 1, and 7 = 3 The solid symbol presents the histogram, using the 
exponential bin plot. 

vanish as the existence of hub is growing larger. In case of small- world network, 
the rewiring probability approaches a threshold to eliminate memory effect, 
the exponent r also follows the mean- field value for different T r . 

The probability distribution of first and all return time is valuable quanti- 
ties for investigating the spatiotemporal correlation and punctuated pattern 
[11,12,13,19]. The definition of first return time with the size t is a separating 
intervals which activated subsequently from a given active neuron. In Fig.5, 
we present first return time distribution(FRTd) for 7 = 3, T r — 1. The first re- 
turn time distribution do not satisfy power-law behavior contrary to the lattice 
case. The power-law behavior of the early return time region is mostly affected 
by the dynamics of a hub node as the hub and nearest neighbors frequently 
fire in the ratio of khub + 1 ,which khub is the degree of the hub and (khub) = 45 
in Fig.5, that is, the hub has a probability of the firing as l/(k hub + 1). An 
interval of power-law regime in early return time increases with the degree of 
the hub. The intermediate return time distribution become almost constant 
for all nodes as iV — > 00 because each node has a same probability to be ac- 
tive again since UCM do not have the degree correlations. Finally, long return 
time has a long exponential decay regime with cut-off for t — > by the finite 
size effect of the dangling node because the diameter of scale-free network is 
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Fig. 6. The log-log plot of the probability distribution function P(s) of the all return 
time as a function of the time t at the critical fitness on UCM with N = 10000, 
m = 3 at T r = 1, and 7 = 3 The solid symbol presents the histogram, using the 
exponential bin plot. 

very small. All return times with the size t is the elapsed time steps to time 
t regardless of the intermediate firing since a given neuron fired at time to- 
rn Fig. 6, we plot all return time distribution for 7 = 3, and T r = 1. The all 
return time distribution(ARTd) also is divided up three characteristic time 
scaling regimes. Furthermore, the slope of each regime is small in early return 
time regime like the lattice model although the general scaling relation is not 
satisfied. 



4 Conclusion 

We have studied the simple model of the neuron function in the brain on the 
uncorrelated scale-free network. Our model show the avalanche dynamics with 
different degree exponent 7 of the networks and different refractory periods T r . 
The avalanche size exponent, r is not changed according to varying 7, for 7 > 
2.4 as well as T r . We measured first and all return time distribution(RTDs). 
The RTDs do not follow power-law behaviors consistently and show the three 
characteristic regions. In future work, the simple model for the SOC system on 
the scale-free network that memory effects contribute would be inquired and 
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a kind of scale-free network more neural system-like should be investigated. 
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